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ABSTRACT 

We determine the magnetic helicity, along with the magnetic energy, at high latitudes using data from the 
Ulysses mission. The data set spans the time period from 1993 to 1996. The basic assumption of the analysis 
is that the solar wind is homogeneous. Because the solar wind speed is high, we follow the approach first 
pioneered by Matthaeus et al. (1982, Phys. Rev. Lett. 48, 1256) by which, under the assumption of spatial 
homogeneity, one can use Fourier transforms of the magnetic field time series to construct one-dimensional 
spectra of the magnetic energy and magnetic helicity under the assumption that the Taylor frozen-in-flow hy- 
pothesis is valid. That is a well-satisfied assumption for the data used in this study. The magnetic helicity 
derives from the skew-symmetric terms of the three-dimensional magnetic correlation tensor, while the sym- 
metric terms of the tensor are used to determine the magnetic energy spectrum. Our results show a sign change 
of magnetic helicity at wavenumber fc « 2 AU^^ (or frequency v ^ 2 /xHz) at distances below 2.8 AU and at 
fc 30 AU^^ (or i/ w 25 /iHz) at larger distances. At small scales the magnetic helicity is positive at northern 
heliographic latitudes and negative at southern latitudes. The positive magnetic helicity at small scales is ar- 
gued to be the result of turbulent diffusion reversing the sign relative to what is seen at small scales at the solar 
surface. Furthermore, the magnetic helicity declines toward solar minimum in 1996. The magnetic helicity flux 
integrated separately over one hemisphere amounts to about 10*^ Mx-^/cyclc at large scales and to a 3 times 
lower value at smaller scales. 
Subject headings: MHD - Turbulence 



I 1. INTRODUCTION 

Over the past 30 years there has been considerable ac- 
tivity in estimating magnetic and current helic ities of the 
Sun's magne t ic fiel d both at t he surface (Seehafed 119901 
■ i Pevtsov st all 119951: iBao et alj [l999; Pevtsov & LatushM 
; 2oooir as well as in the s olar wind (Matthaeus et al., il982t 
, iRust & Kumaiill994l 1199 6) using a variety of spacecraft, in 
particular Voyager 2. The early motivation of lSeehafeJ (11990 1) 
was the connection with the a effect in mean-field dynamo 
theory. Subsequent work confirmed his early findings that the 
current helicity has a negative sign in the northern hemisphere 
and a positive in the southern. This also agreed with expecta- 
tions according to which current helicity is a proxy for kinetic 
helicity (Keinigs 1983), which is known to be negative for cy- 
clonic events in the northern hemisphere and positive in the 
southern. 

Later work in connection with dynamo theory in periodic 
domains clarified that a correspondence between kinetic and 
current helicities can only be expected to hold at and below 
the scale of the energy-carrying eddies of the turbulence, be- 
cause at larger scales the signs of current and magnetic helic- 
ities should reverse ([Brandenburg 2001). This is connected 
with magnetic helicity evolution and the fact that in a effect 
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dynamos magnetic h elicity at large a nd small scales tend to 
have opposite signs (lSeehafejll996t ll [T999h . If magnetic 
helicity fluxes and resistive effects are weak or unimportant, 
e.g., in the kinematic regime of a growing dynamo, the to- 
tal magnetic helicity is constant or zero if it was zero ini- 
tially. The production of magnetic helicity at the scale of the 
energy-carrying eddies is then accompanied by the production 
of magnetic and current helicity of the opposite sign at scales 
larger than the scale of the energy-carrying eddies. There is a 
similar tendency also when total magnetic helicity is not con- 
served, e.g., on long timescales when resistive effects become 
important or if magnetic helicity fluxes are present. Thus, 
there are good theoretical reasons to expect that the magnetic 
field in the Sun is bi- helical, i.e., of opposite sign at s mall and 
large length scales ( Blackman & Brandenburg|2003h . 

We emphasize that 'small' refers here to the scale of the 
energy-carrying eddies, which is also called the outer scale — 
in contrast to the inner scale where kinetic and magnetic ener- 
gies get dissipated. In the Sun the outer scale can be as large as 
50 Mm, which corresponds to the pressure scale height near 
the bottom of the solar convection zone. On the other hand, 
'large' refers to the scale of the mean field that characterizes 
the solar cycle. The width of the toroidal flux belts, i.e., the 
width of the wings of the butterflies in a solar butterfly di- 
agram is about 20°, corresponding to about 200 Mm in the 
Sun. In any case, the scale of the large-scale field is finite, i.e., 
such a field is still subject to decay and possible regeneration 
by dynamo action, and should thus not be confused with an 
'imposed' magnetic field. The latter case is sometimes con- 
sidered in numerical simulations, where the departure from an 
imposed field corresponds to the small-scale field whose heli- 
city i s not conserved on its own ('Stribli ng et alJ[T995t iBergeii 
119971; [Brandenburg & Matthaeu s 2004). Let us also mention 
at this point that the magnetic helicity is quadratic in the mag- 
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netic field, so it is not expected to flip sign from one cycle to 
the next, although it may of course vary in strength. 

The bi-helical nature of the ma gnetic field has been the 
topic of related work by I Yousef & Brandenburg ( 2003 ), who 
investigated the relaxation of an initially bi-helical field 
and the mutual annihilation of the two signs of magnetic 
helicity. It should be noted that a connection has also 
been discussed between the current helicity observed in 
the Sun and that obtained from mean-field dynamo mod- 
els (Seehafer 1^9^, "M^ati&iGnmiBllOOl!). Furthermore, 
IChoudhuri et al. (2004) find mostly negative current helicity 
in the north, except that during short intervals at the begin- 
ning of eac h cycle the cu r rent h elicity in the north can be pos- 
itive, while IZhang et"an (120061) find a band of negative cur- 
rent helicity at mid-latitudes and positive values at higher and 
lower latitudes. However, these papers ignored the possibil- 
ity that the magnetic field in each hemisphere is expected to 
be bi-helical. The first observational evidence for a reversed 
sign of magnetic helicity at large scales came from an anal- 
ysis of synoptic maps of the radial magnetic field of the Sun 
dBrandenburg et alj2003h . They found a sign reversal of mag- 
netic helicity at the time of solar maximum with positive val- 
ues after that moment. 

There is yet another reason for studying magnetic helicity 
in the Sun. Magnetic helicity is a topological invariant that is 
equal to half the number of constructive flux rope c rossings 
times the square of the magnetic flux in these ropes (" Moffatti 
[1969). Magnetic helicity is therefore a measure of the de- 
gree of tangledness. It might then be possible to assess the 
degree of tangledness by countin g the net cr ossings of fila- 
ments in Hq images of the Sun (iChaell2000h . although this 
technique leaves some ambiguity regarding the sign of the 
magnetic helicity. Other ways of measuring magnetic helicity 
and their fluxes is by tracking the motions at the solar sur- 
face ( Kusano et al. 2002; Demoulin & Berger 2003), which 
led to the estimate that the total flux of magnetic helicity in 
each hemisphere, integrated over a full 1 1 -year cycle, is of the 
order of 10''*' Mx^/cycle dBerger & RuzmaikinlllOOOl) . This 
number agrees also with theoretical expectations of an upper 
limit of this value (Brandenburg & Sandin 2004; Brandenburg 
[2009). However, there is no indication as to a possible scale 
dependence of the magnetic helicity. The only evidence for 
this is just the qualitative appearance of a systematic tilt of 
bipolar regions. This tilt corresponds to writhe helicity, which 
is a quantity that depends only on the topology of the axis 
of a flux tube structure. Independent of the time during the 
1 1 -year cycle, it should have a positive sign in the northern 
hemisphere and a negative sign in the southern hemisphere, 
i.e., just the opposite of what is observed in the magnetic field 
line twist at smaller scales. 

The hope is now that measurements of the solar wind might 
help teach us something about the scale dependence of the 
contribution to the magnetic field that is related to the a ef- 
fect. In order for the a effect to work efficiently and to 
escape what is known as catastrophic a quenching, nega- 
tive magnetic helicity associated with small-scale m agnetic 
fields must be shed (B lackman & Field 2000; Kleeori n et al.l 
|2000() ; see [Brandenburg & Su bramanian (2005) for a re- 
view. This might be accomplished b y coronal mass ejec- 
tions dBlackman & Brandenburg|l2003l) . Coronal mass ejec- 
tion events are manifol d, and they are almost all associated 
with magnetic helicity (iDemoulin et al.l [20021) . but concern 
only the corona and not the solar wind. The large-scale mag- 



neti c field in the solar wind is characterized by the Parker spi- 
ral (lParkerllI958l) . The helicity associated with the Parker spi- 
ral is known to be negative in the northern h emisphere and 
positive in the southern (Bieb eret al.l II 98731) - independent 
of the time during the 11 -year cycle. Therefore, even though 
we would normally associate the Parker spiral with the large- 
scale field, its helicity is of opposite sign to the helicity of the 
large-scale field generated in the dynamo interior, or that ex- 
pected from the tilt of the flux tubes near the solar surface. 
On the other hand, the sign of the helicity associated with the 
small-scale field that needs to be shed does agree with that of 
the Parker spiral, although it would seem to be of the wrong 
scale. Nevertheless, not much is known about the relationship 
between magnetic helicity fluxes and magnetic helicity itself. 
For example, it is possible that turbulence in the region out- 
side the dynamo would continue to diffuse the magnetic field, 
although it would no longer amplify it by an a effect. This 
effect would tend to reverse the production of bi-helical mag- 
netic fields and would pump positive magnetic helicity into 
smaller scales, leaving behind negative magnetic helicity at 
larger scales. This could be interpreted as a forward turbulent 
cascade, but it is probably only possible in an expanding flow, 
so as to not cause a conflict with the realizability condition 
that enforces the inverse transfer in a confined helical flow. 
This phenomenon was seen in me an-field calculations with 
a turbulent exterior (see Fig. 7 of [Brandenburg et al J 120091) 
and also in direct numerical simulations of dynamos in spher- 
ical geometry with a n early force-free exterior (see Fig. 4 of 
iWarnecke et"al]|20I lal) . 

To probe quantitatively the possible scale dependence of the 
magnetic field, perhaps the best type of analysis is that used in 
the early meas urements on Voyager 2. Making use of the Tay- 
lor hvpothesis. lMatthaeus et al.] (II982I) were able to associate 
frequencies with wavevectors. Making the further assumption 
of homogeneity, they were able to translate the simultaneous 
measurement of the two field components perpendicular to 
the direction of the wind to information not only about the 
magnetic energy spectrum, but in particular, about the mag- 
netic helicity spectrum. The background for application of 
this technique to the computation of other helicities was ex- 
plored further in Matthaeu s et al.l (11986a). 

A possible complication with many of the early results is 
that the trajectories of Yoyager and other spacecraft were close 
to the ecliptic, across which the magnetic helicity is expected 
to change sign. However, there are no published spectra that 
appear to involve data sets that crossed the heliospheric cur- 
rent sheet. Nevertheless, the magnetic helicity as found by 
Matthaeus & Goldstein ( 1982) randomly changed sign at all 
scales, although Goldstein et al] dI99Ih did find short inter- 
vals during which the magnetic helicity had a constant sign 
at scales close to that of the proton gyroradius. Further- 
more, Smith & Bieber (1993) found that at frequencies be- 
low 10 /iHz the magnetic helicity tends to have a predominant 
sign: negative in the north and positive in the south, while at 
higher frequencies there are strong fluctuations of the sign. 
However, we shall argue below that, by averaging over broad 
wavenumber bins, it is still possible to extract meaningful in- 
formation from the data even at high frequencies. The tech- 
nique applied by Matthaeus et al. ( 1982) appears very suitable 
for the purpose of assessing scale dependence of magnetic he- 
licity. The purpose of the present work is therefore to apply 
this technique to more recent measurements of Ulysses that 
flew in a nearly polar orbit that covered both hemispheres. 
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2. DATA ANALYSIS 

We use 60 s time averages from the Vector Helium Mag- 
netometer on Ulysses. The original time resolution is up 
to 2 vectors/second and the sensitivity is about 10 pT; see 
iBalogh et alj (119921) for a detailed description. The available 
data comprise measurements of all three components of the 
magnetic field B and velocity u in the locally Cartesian helio- 
spheric coordinate system {R, T, N), where R is the distance 
from the Sun, T points in the transverse direction parallel to 
the solar equatorial plane and is positive in the direction of 
solar rotation, and N = Rx T is the third direction pointing 
toward heliographic north. This corresponds to a right-handed 
coordinate system. Note that the R-T plane is inclined to the 
heliographic equatorial plane by an angle equal to the helio- 
graphic latitude A of the spacecraft. 

We have analyzed 27 data sets comprising a time span of 
about one month each and covering different epochs between 
1993 and 1996. During 1993/94, Ulysses was at 39° to 80° 
southern latitudes and distances between 4.3 AU to 1.7 AU, 
while during 1995/96 it was at 43° to 79° northern latitudes 
and distances between 1 . 5 AU to 3 . 5 AU. For each data set we 
determined the average radial wind speed un, which ranges 
between 720 and 790km/s. Note that the local wind speed 
points almost exactly in the direction away from the Sun. Us- 
ing Taylor's hypothesis, we translate time t into the negative 
radial coordinate, R ^ Ro — upt, where i?o is the slowly 
changing distance of the spacecraft. Next, we compute the 
Fourier transform of each of the field components. 



^""B.iR) dR, R,T, N, (1) 



and are thus able to compute the spectral correlation matrix. 



M"^(fcfl,)=i?,(fcfl,)B*(fcfl) 



(2) 



where the asterisk denotes complex conjugation. The su- 
perscript ID emphasizes an important difference with the 
three-dimensional correlation tensor Mfp{k). Before dis- 
cussing this in more detail we note that, in practice, AllJ^ikfs) 
is obtained from measurements along the R direction. In 
that case one computes the one-dimensional magnetic en- 
ergy and helicity spectra simply as jii^E^ {ka) = \IS\'^ and 

HhP(kR) = ilmiBrBtr ) /fc/?. These are the equations used 
by lMatthaeus et al](ll982h . and we shall use them for the most 
part of this paper as well. 

It is well known that no assumption about isotropy is made 
in obtaining the one-dimensional magnetic energy and heli- 
city spectra. However, we should emphasize that, even if 
the turbulence were isotropic, the one-dimensional spectra 
obtained through direct measure ments are not equival ent to 
the three-dimensional ones (T ennekes & Lum ley 1972). The 
differences can become important in regi ons where the spec - 
tra deviate from pure power law scaling jDobler et al.ll2003h . 
The purpose of the rest of this section is to extend the well- 
known formula for the conversion between one- and three- 
dimensional energy spectra to the case of helicity spectra. To 
highlight the analogy between the two, we ignore here the 
consideration of longitudinal and transverse energy spectra 
and make the assumption of isotropy. This assumption does 
seem at odds with results obtained in the ecliptic (see, e.g., 
iNarita et al.l2010t[S^aoui et al.l20 10'), but is consistent with 
an analysis of Ulysses data reported by .Smith (.2003.) . Even 



though there is near isotropy of the variances (shown also be- 
low), there is no spectral isotropy, so the correlation length 
perpendicular to B is shorter than along B. However, by 
making the assumption of isotropy, we shall be able to assess 
the differences between three- and one-dimensional spectra. 
It will turn out that these differences are rather small. 
In three-dimensional isotropic helical turbulence we have 



2fioEl^{k) ikkH^ik) 
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where k — k/kis the unit vector of k. Note that these spectra 
obey the realizability condition. 



2f,oEl^ik)>k\HiP{k)l 



(4) 



where the factor 2 in front of Em {k) is just a consequence of 
the factor 1/2 in the definition of energy. The two spectra are 
normalized such that 

J 5,jMf{k) d^A: = f^oE^ik) dk = {B^)/2, (5) 

/•I nOO 
^-jik^MfJ'ik) d^k = HiPik) dfc = (A . B), (6) 

where JB = V x A is the magnetic field expressed in terms 
of the magnetic vector potential A, which obeys the Coulomb 
gauge, V • A = 0, and angular brackets denote averaging over 
the data spanned by each data set. 

Let us now relate AffP(fc) to M/P(fci?). Suppose MfJ^{k) 
were known, then, to improve the statistics, one can obtain 
Mlj^{kfj) from MfJ^{k) by averaging over the other two 
wavevector components, i.e.. 



MlJ'ikR) 



Mff{kR, kT,kN) dkr dk^- 



(7) 



We define one-dimensional energy and helicity spectra via 



5,,MlJ'{kR)^^jiaE^{kii), 



mi 



2e^jR{xkR)-HdlJ'{kR) = Hl^{kR) 



ID/ 



(8) 
(9) 



and note that they are related to the three-dimensional energy 
and helicity spectra via an integral transformation 



El]^{kR) 



EiY{k)d\-ak, 



HlJ^ikR) = / i/^^(fc)dlnfc 



(10) 



(11) 



This transformation is well known for the e nergy spectrum 
(cf., iTennekes & Lumlevl[T97l iDobler et al.l 2003), but has 
been generalized here to the case with helicity; see Ap- 
pendix |A] for details of the derivation. In the following we 
present results first for E^ and and compute then the 
three-dimensional spectra via differentiation, i.e.. 



E^{k) = -dEl^{k)/d\nk, 



-dHl^{k)/d\nk. 



(12) 
(13) 



Note, however, that differentiation amplifies the error in the 
already rather noisy data. Therefore we perform differentia- 
tion based on data that have been averaged into rather broad 
wavenumber bins. We use a second-order midpoint formula 
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with respect to In A: bins. This yields data at k values that 
lie between those for the one-dimensional spectra. In the fol- 
lowing, where we present mostly one-dimensional spectra, we 
omit the superscript ID, but retain superscript 3D for all three- 
dimensional spectra. 

3. RESULTS 

For the 27 data sets analyzed, the distance of the space- 
craft to the Sun varies from ^ 1.5 AU to ^ 4.5 AU. This 
needs to be taken into account when combining different data 
sets. In Figure [T] we show that the integrated magnetic en- 
ergy in the three field components is approximately equal 
in the parameter range covered by the 27 data sets. Th is is 
compatible with earlier results of lMatthaeus et al.l (Il986bl) us- 
ing data from Voyager 2. Furthermore, all three contribu- 
tions fall off slightly faster with distance than oc R^^. One 
would expect a perfect Br^ scaling if the Poynting flux stayed 
constant, which one might expect for a magnetically domi- 
nated wind. This suggests that magnetic energy is dissipated 
into heat, which has been discussed in detail in recent years 
fcoldstein et al. 1995; Tu & Marsch 2003; Freeman 1991 
[Smith et al. 2001; S ahraoui et al.ll2009L 12010 ). An estimate 
for the corresponding magnetic 'luminosity', assuming ap- 
proximate isotropy over the solid angle, i.e.. 



u-dS ^47rR^{B^/2na)uR, (14) 



also falls off from 2.8 x IQi^ W at 1.5 AU to 1.2 x lO^^ W 
at 4.5 AU. The magnetic luminosity based on the perpen- 
dicular field components shows a weaker decline from 1.6 to 
0.8 X 10^^ W. This might be a consequence of th e Parker spi- 
ral for which one expects a Bt / Br oc R scahng (iBieber et al.l 
ll987atfWebb et al.l20Toh . The present data suggest that this is 
however a weak effect. Our proxy for the magnetic luminosity 
corresponds to (2-4) x IQ^^Lq, where Lq is the bolometric 
luminosity of the Sun. Thus, because of the approximate R^^ 
dependence we scale the spectra Eyiik) and Huik) to a ref- 
erence distance of 1 AU before averaging over different data 
sets. 

Next, we focus on the magnetic energy and helicity spec- 
tra. We begin by presenting results where we combine data 
from both hemispheres. As we will verify later, the magnetic 
helicity has opposite signs in the northern and southern hemi- 
spheres, so we multiply the helicity measured in the south 
by —1. Another possibility would be to divide by sin A, but 
at least for the Parker spiral one does expect a much sharper 
sign change near the equator tha n what is expected from a 
sin A profile (iBieber et al.lll98'7ah . Furthermore, we distin- 
guish between data sets where the distance to the Sun is ei- 
ther inside or outside 2.8 AU. In Figure |2] we plot, separately 
for two separate distance intervals, 2/ioi?M(fc) and k\HM{k) \, 
rescaled by AttR^, as well as the relative magnetic helicity, 
kHM{k)/2iioEM{k). We also show a cumulative average of 
this ratio, starting from the low wavenumber end. This shows 
quite clearly that the magnetic helicity in the north is neg- 
ative at small wavenumbers (large length scales) and that it 
becomes positive at large wavenumbers (small length scales). 
The cumulative nature of the average slightly overemphasizes 
the negative contributions, delaying thus the position where 
the sign changes. The negative magnetic helicity at large 
scales agrees with that belonging to the Parker spiral, while 
the positive magnetic helicity at small scales could be the re- 
sult of turbulent diffusion leading to a reverse transfer of mag- 
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Fig. 1 . — Magnetic energy density and magnetic luminosity as a function 
of distance. Note that the decay of magnetic energy is slightly faster than 
oc and that the magnetic 'luminosity' varies in the range 4-8 X 10^^ W. 

netic helicity from large to small scales, as explained in the 
introduction. 

In the remainder of this paper we use logarithmically 
spaced wavenumber bins. This way we can reduce the data 
to a small number of bins and thereby minimize the statisti- 
cal noise in a more meaningful way. This is shown in Fig- 
ure [3] In each bin we average the actual values (not the 
logarithms) of spectral energy and magnetic helicity, both 
weighted with a k factor. We use rather broad bins, for ex- 

T-l 



1.5 AU has a width 

;AU"^ has 



ample the first wavenumber bin at k 
of Ak — 1.2 AU^^, and the second bin at = 5. 
Ak = 4.8AU"\ Given a maximum length of 1 month for 
each of the 27 data sets, we have in principle a spectral reso- 
lution of ACl = 0.2d"\ corresponding to Afc = 0.5AU"\ 
We use data that were already averaged over 60s time in- 
tervals. This corresponds to a spatial resolution of 50 Mm 
and hence a Nyquist wavenumber of 10"* AU^^. At distances 
below 2.8 AU, the magnetic helicity is negative only in the 
smallest wavenumber bin (k « 1 AU^^), while at distances 
beyond 2.8 AU the first 3 wavenumber bins (fc < 30 AU"^) 
show negative helicity. 

An estimate for the error has been obtained by compar- 
ing with the averages that result by taking only data from the 
northern or the southern hemisphere into account. The larger 
one of the two departures is taken as an estimate of the error 
This results in a relative uncertainty of our average values by 
a factor of 1-2. 

As explained above, three-dimensional spectra of magnetic 
energy and magnetic helicity can be obtained by differenti- 
ation with respect to Ink. It turns out that, within expected 
error margins, the three-dimensional spectra are surprisingly 
close to the one-dimensional spectra. This is shown in Fig- 
ured where we compare the types of spectra. In view of the 
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Fig. 2. — Magnetic energy and helicity spectra, 2/^o^M(fc) and 
k\HM{k)\, respectively, for two separate distance intervals (first and third 
panels). Furthermore, both spectra are scaled by AnR^ before averaging 
within each distance interval above and below 2.8 AU, respectively. The rel- 
ative magnetic helicity, kH-^{k) /2^QE-!^i(k), is plotted separately (second 
and fourth panels) together with its cumulative average starting from the low 
wavenumber end. The zero line is shown in dashed. 



fact that the three-dimensional spectra do not seem to alter 
our conclusions, and since differentiation increases the noise 
in the data, we restrict ourselves in the following to the dis- 
cussion of one-dimensional spectra. 

It is worthwhile noting the large separation between both 
graphs in the ordinate. In other words, k\HM{k)\ <^ 
2/ioi?M(fc)- This indicates that the relative magnetic helicity 
is rather small, which may not be too surprising considering 
the fact that we have averaged a noisy magnetic helicity result 
over rather broad wavenumber bins. Also, of course, there is 
no reason to expect the relative magnetic helicity at the so- 
lar surface to be particularly strong. Next, from the open and 
filled symbols we can see the sign of Huik), which turns out 
to be negative at the largest scales (filled symbols) and posi- 
tive at small scales (open symbols). At a larger distance from 
the Sun, the break point where the sign of the magnetic heli- 
city changes grows to larger wavenumbers, corresponding to 
smaller scales. This is probably again related to the effect of 
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Fig. 3. — Magnetic energy and helicity spectra, 2fioEiji{k) and kHyi{k), 
respectively, for two separate distance intervals. Fuilhermore, both spectra 
are scaled by AttR^ before averaging within each distance interval above and 
below 2.8 AU, respectively. Filled and open symbols denote negative and 
positive values of Hyi{k), respectively. 
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Fig. 4. — Comparison between 3D and ID spectra of magnetic energy and 
helicity. In both panels the ID spectra are denoted by dotted lines, while the 
3D spectra of magnetic energy by dashed and the 3D spectra of magnetic 
helicity are indicated by filled and open symbols for negative and positive 
contributions. 

turbulent diffusion causing the reversed transfer of magnetic 
helicity from larger to progressively smaller scales. It could 
also be a consequence of the growing dominance of the large- 
scale field (having negative helicity) with distance, so that it 
would appear as if the magnetic hehcity of the large-scale 
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Fig. 5. — Magnetic helicity (open symbols) and energy (filled symbols) vs. 
latitude for data from wavenumber bin k = 1.2 AU~^ and distances R < 
2.8 AU (upper panel) and wavenumber bin k = 300 AU~^ and all distances 
(lower panel). In the last panel the energy is downscaled by factor 10, because 
otherwise those data points would lie outside the plot range. The solid lines 
represent fits proportional to sin A, while dash-dotted lines represent the zero 
value. 

field imprints itself onto the smaller scales. A similar effect 
has been seen in helical dy namo simulations with an irnpose d 
magnetic field (see Fig. 3 of lBrandenburg & Matthaeusl2004l) . 
where for sufficiently strong fields the sign of the magnetic 
helicity is equal to that of the kinetic helicity at small scales. 
At the same time, the relative magnetic helicity diminishes, 
which is indeed also seen in Figure [3] 

In view of dynamo theory and for comparison with earlier 
work, it is of interest to compute magnetic helicity fluxes sep- 
arately for large and small length scales, and integrate them 
over half the solid angle and over the 1 1 year cycle, Tcyc, i.e., 
we compute 



2 -^H eye 



eye 



HM{k)dk, (15) 



where the 1/2 factor takes the fact into account that one 
normally gives magnet ic hehcity fluxes integrated separately 
for each hemisphere (iBerger & RuzmaikinI I200Q) . In Ta- 
ble [T] we give the results for ^C^T^yc separately for large 
(— ) and small (+) scales and also for small and large dis- 
tances. It turns out that these values are typically around 
10^^ Mx^/cycle, which is remarkably close to early estimates 
of lBieber & Rus3 (119951) of 2 x lO^^ Mx^ /cycle, and about 10 
times below the expected upper limit (iBrandenburg.2009i) . 

TABLE 1 

Results for ^C^Tcyc in units of Mx^ /cycle. 



Distance 



Large Scales 



Small Scales 



B. < 2.8 AU 
R > 2.8 AU 



-0.9 X lO*'' 
-1.3 X 10*^ 



35" 



+0.3 X 10 
+0.03 X 10*5 



Fig. 6. — Similar to Figure|5] but vs. time. Again, in the last panel the 
energy is downscaled by factor 10, because otherwise those data points would 
lie outside the plot range. The dash-dotted lines represent the zero value. 

Next, we consider the latitudinal dependence by abandon- 
ing the averaging over heliographic latitude and consider 
data from two separate wavenumber bands around 1.2 and 
300AU~\ The data are obviously very noisy now, espe- 
cially at low wavenumbers where the wavenumber bins in- 
volve fewer data points; see Figure |5] Nevertheless, there is 
still some evidence for the magnetic helicity having opposite 
signs in the two hemispheres and, in addition, opposite sign at 
large (5 AU) and smafl (0.02 AU « 3000 Mm) length scales 
(= 27r/fc). 

The present data have all been taken from the time just af- 
ter solar maximum and before the next solar minimum. It 
would therefore be interesting to compare this with measure- 
ments taken at other times. There is in principle the possibility 
that the magnetic helicity could change and even reverse sign 
for brief time intervals. To see whether already the present 
data suggest a possible systematic temporal trend, we plot the 
spectral magnetic helicity in the same two wavenumber bins 
(around 1.2 and 300 AU^^) vs. time; see Figure|6l Given the 
high noise level, it is not possible to conclude anything def- 
inite from the plot. Only at large wavenumbers (shown here 
for k = 300 AU'^) there might be a meaningful downward 
trend. Of course, one can only be really sure about this when 
comparing data over a much longer time span, in which case 
one might expect to see an oscillatory variation as is also seen 
in Figures 3 and 4 of.Ber ger & R uzmaikin (2000,) . 

4. CONCLUSIONS 

The present results have revealed for the first time evi- 
dence that the magnetic helicity in the solar wind has op- 
posite signs at large and small length scales. However, 
the signs are actually the other way around than what was 
naively predicted based on th e expected signs at the solar 
surface dBlackman & Branden burg 200 3!), but they d o agree 
with more recent simulations of Brandenburg et alj ( 12009b 
and Warnecke et al.. (2011a.b.) that show a reversed sign some 
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distance away from the dynamo regime due to the effect of 
turbulent diffusion (without a effect) that tends to forward- 
cascade magnetic helicity from large to small scales. 

Although the solar wind would carry some imprint of 
the fields generated within the Sun, this field would be de- 
formed due to the fact that the material leaves the Sun 
with a net angular momentum, and is perhaps turbulent. 
The magnetic helicity generation by turbulence may be un- 
derstood using the magnetic helicity conservation equation 
for the gauge invariant small-scale helicity h[ derived by 
Uubr amanian & Brandenburg (2006). We have akeady used 
this argument in this paper when we associated the effect of 
turbulent magnetic diffusion with a forward turbulent cas- 
cade. Suppose that the wind is turbulent on small scales, 
then the helicity generation on those scales is governed by 
dhf/dt — —2£ ■ B + where £ is the mean turbulent 
electromotive force and the dots refer to microscopic dissipa- 
tion and p ossible flux terms. As in mean-field d ynamo theory 
(see, e.g.. [Brandenburg & Subramanianll2005h . we approxi- 
mate £ ^ aB — rjtJ, where J = V x B is the mean 
current density, and a and % refer to a possible a effect 
and turbulent diffusion associated with solar wind turbulence, 
then dhf/dt — —2aB^ + 2ritB ■ J. The dynamo gener- 
ated large-scale field B, which has positive helicity, will then 
lead through the rjt term above to small-scale helical fields 
also with positive h[. Moreover, because total helicity is con- 
served, the turbulent diffusion will put a negative helicity h^^ 
on larger scales. Since this ft-m was originally positive, this 
leads to a decrease of h^. This process may be behind the 
apparent forward cascade of helicity. 

Now consider the large-scale field in the solar wind. Re- 
alistically, this field will also have reversals associated with 
the poloidal field from the Sun reversing every 11 years. 
However, a unique sense (sign) of its helicity can be im- 
printed because the sense of rotation of the Sun is always 
the same, and, in addition, the radial velocity is a l ways ou t- 
ward. We know from the work of Biebe r et al.l ( Il987al[bl) 
that this sense is to give negative helicity in the north, op- 
posite to what the dynamo generated mean field had in the 
solar interior and at the surface. Perhaps the differential rota- 
tion in the solar wind pumps negative helicity to the north 
and positive to the south so as to reverse the original sign 
of the mean field helicity. One can then explain our current 
observations as follows: First, the positive helicity seen at 
larger k is that generated from the turbulent diffusion spread- 
ing the helicity of the dynamo-generated mean field to larger 
k and then advecting this outwards. Second, the Parker spi- 
ral eventually leads to a negative helicity on the largest scales 
dBieber et al.|[T98 7a). In this picture, the helicity of the Parker 
field would correspond to scales around 5 AU (corresponding 
to fc = 1.2 AU"^ at i? < 2.8 A U; see Figure |3]l. This is also 
in agreement with early work of 'S mith & BiebeJ (11993 ) who 
found negative magnetic helicity in the north below frequen- 
cies of about 10 ptHz, corresponding to fc < 20 AU'^ at a 
wind speed of about 400 km/s. However, this scenario needs 



to be much better explored through simulations such as those 
of Wa rnecke et al.1 (1201 IH) . 

Unfortunately, the relative magnetic helicity is rather weak. 
Therefore, only through extensive averaging we are able to 
extract any useful information. This low level of relative 
magnetic helicity suggests that there is efficient mixing tak- 
ing place in the solar wind, but it could also mean that the 
magnetic helicity is akeady rather low at the solar surface. 
Another possible reason could be the proximity to the solar 
minimum in 1996. 

Although magnetic helicity fluxes are expected to keep their 
preferred sign over the solar cycle, some mo dulation is defi- 
nitely to be expected, as was already found bv Smith & Bieberl 
( 119931) using measurements spanning the years 1965 to 1988. 
The data span analyzed in the present work is too short to 
make any meaningful statements, but one can see that at least 
at larger wavenumbers there seems to be a decline in mag- 
netic helicity as one progresses further toward solar minimum 
in 1996. 

The present data can be used to extract quantitative esti- 
mates for magnetic helicity fluxes. This quantity is normally 
quoted in Maxwell squared per solar cycle. Assuming that 
the magnetic helicity during 1993-1996 is representative of 
the rest of the solar cycle, our analysis suggests values around 
10*^^ Mx^ /cycle, wh ich is comparable to the earlier work of 
iBieber et al.1 ( I1987H) . 

Although the results presented in this paper are physi- 
cally appealing, there remains uncertainty about the assump- 
tions made in this work, most notably the Taylor hypothe- 
sis, which may be problematic over large length scales, but 
it is expected to be reasonably well satisfied for fluctua- 
tions in the inertial range of the turbulent spectrum (see, e.g., 
lMatthaeuF& Goldstein 1982). The assumption of isotropy 
is only needed when computing three-dimensional spectra. 
Only the two components in the plane perpendicular to the 
radial direction enter in our analysis of magnetic helicity. The 
finite pitch angl e of the Parker sp iral introduces anisotropy in 
that plane (Matt haeus et al.|[T996l) . T his becomes important at 
progressively smaller length scales (IWicks et al.ll2bl ll) . even 
though the variance of the three components remains similar 
(see Figure [T)- Using spe ctra of the different components of 
the magnetic field, Bi eber et al.l (Il996i) was able to determine 
that 85% of the magnetic energy resi des in the perp endic- 
ular field fluctuations, while for Ulysses, ISmithI (120031) found 
that this value is closer to 50%. In any case, it would be use- 
ful to estimate magnetic helicity using synthetic data from a 
numerical simulation. Such an exercise is likely to provide 
useful insight into the reliability of the assumptions made and 
the accuracy of the method. 
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APPENDIX 

DERIVATION OF EQUATIONS (10) AND (II) 

We begin with Eq. (|7]i, but instead of performing the integration over fc^ and k^ we use cyUndrical coordinates {k±, (pk) in 
Fourier space and assume axisymmetry, so therefore dkrdkN = 2TTk±dk±. Thus, we have 

MlJ'ikR)^ / MfJ'{kR,k±,cPk)k±dk±dcj)k^2TT M!°{kR,k^)k±dk±. (Al) 
Jo Jo Jo 
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Next, we insert Eq. (O, take the trace, use — kj^ + k\, and substitute k^dk±^ — kdk valid for a fixed kji, to carry out the 
integration over k in the allowed range from ku to cxd, and obtain, using Eq. ([8]), 

6,,M!^{kn) = 2n P ^ g W ^ ^fe = H ^"^mW ^ H ^,El^(k)dlnk, (A2) 
which corresponds to Eq. ( fTOl ). Likewise, multiplying instead with —2eij jiikj^^ , and using Eq. (|9]l, we arrive at 

2e,,n{^kn)-'Mip{kn)^27r H A^^^^kdk^ H EER ik = H H!Eik)d\nk, (A3) 

which corresponds to Eq. (fTTT i. In particular, this yields fJ,oEl^{kii) = \B\'^ and Hl['{kii) — 4Im(i?TS^)/fcfl, as stated just 
below Eq. 
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